Dicke-like effect in spin-polarized transport through coupled quantum dots 
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Spin-dependent electronic transport through a quantum dot side-coupled to two quantum dots and 
attached to ferromagnetic leads with coUinear (parallel and antiparallel) magnetizations is analyzed 
theoretically. The intra-dot Coulomb correlations are taken into account, whereas the inter-dot ones 
are neglected. Transport characteristics, i.e. conductance and tunnel magnetoresistance associated 
with the magnetization rotation from parallel to antiparallel configurations, are calculated by the 
noneqiulibrium Green function technique. The Green functions are derived by the equation of 
motion method in the Hartree-Fock approximation. The conductance spectra are shown to reveal 
features similar to the Dicke resonance in atomic physics. 
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I. INTRODUCTION 

Interference effects in electronic transport through cou- 
pled quantum dots (QDs) are of great interest from both 
fundamental and application points of view. One of the 
simplest systems where the interference effects are clearly 
visible, consists of two QDs connected in parallel to elec- 
tron reservoirs (metallic or semiconducting leads). This 
interferometer-like geometry allows observation of such 
quantum phenomena like the Aharonov-Bohm oscilla- 
tions and the Fano antiresonance^'^ii^iii^"ii^"Sii2,iliiiSiiil. 

Conductance of systems consisting of three or 
more QDs reveals further interesting features of 
quantum transporO^^iii^iia^^^^^. Of par- 
ticular importance seems to be the interference 
phenomenon2^i^i22ii^i2^, which resembles the well known 
Dicke resonance in atomic physics^i^i^. The key fea- 
ture of the Dicke resonance is the presence of a strong and 
very narrow emission line of a collection of atoms which 
are separated by a distance smaller than the wavelength 
of the emitted light^Si^. 

In this paper we consider the Dicke-like resonance in 
spin polarized transport through a system of three cou- 
pled QDs. One of the QDs is attached to two ferromag- 
netic leads, while the other two dots are side-coupled to 
the first one. Our considerations include the intra-dot 
electron correlation, while the inter-dot Coulomb inter- 
action is assumed to be much smaller than the intra-dot 
one and is omitted. Transport characteristics, in partic- 
ular conductance of the system and tunnel magnetore- 
sistance associated with rotation of the electrodes' mag- 
netizations from antiparallel to parallel alignment, are 
calculated by the non-equilibrium Green function tech- 
nique in the linear response regime. Since the systems 
with Coulomb interactions cannot be solved exactly, we 
applied the Hartree-Fock (HF) approximation scheme to 
calculate the relevant Green functions from the equations 
of motion. 

In section 2 we describe the model, while theoretical 
formulation of the problem and basic analytical formula 
are presented in section 3. Numerical results for symmet- 



ric and asymmetric systems are shown and discussed in 
section 4. It is shown there that narrow Dicke-like peaks 
appear in the linear conductance. Final conclusions are 
given in section 5. 



II. DESCRIPTION OF THE MODEL 

In this paper we consider a system consisting of three 
single-level quantum dots, coupled as shown schemati- 
cally in figure 1. The dots QDl and QDS are coupled 
to the dot QD2 via direct hopping term. The dot QD2, 
in turn, is additionally attached to ferromagnetic leads 
(see figure 1). Generally, couplings between the dots 
as well as those between the dots and external leads 
can be controlled by applied gate voltages. We also 
note that technological difficulties associated with at- 
taching a semiconductor quantum dot to ferromagnetic 
metallic leads have been overcome and first experimen- 
tal results on spin polarized transport through individual 
semiconductor quantum dots have been published very 
recently'^ . Thus, experimental investigation of spin po- 
larized transport in systems of coupled quantum dots 
seems to be only a matter of time. For simplicity we 
consider only coUinear (parallel and antiparallel) mag- 
netic orientations of the leads' magnetic moments. Apart 
from this, we include the intra-dot Coulomb interactions, 
while the inter-dot Coulomb repulsion is assumed to be 
weak and therefore is neglected. The system under con- 
sideration is then described by Hamiltonian of the gen- 
eral form H = H]^ + H^i + Hqd + Ht- The terms 
Hi^ and Hn describe respectively the left and right elec- 
trodes in the non-interacting quasi-particle approxima- 
tion. Ha = Z^kcr SkacrCkacrCkQfT (for tt = L, R). Here, 



(ckacr) is the creation (annihilation) operator of an 
electron with the wave vector k and spin a in the lead a, 
whereas ekao- denotes the corresponding single-particle 
energy. 

The third term of the Hamiltonian H describes the 
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III. ANALYTICAL DESCRIPTION 

Electric current J flowing through the system is given 
by the formulck^i^i 

+ [fL{e)T\e) ~ fn{e)T^{e)Wie) - G^{e)]}, (3) 



FIG. 1: Schematic picture of the quantum dot system cou- 
pled to ferromagnetic leads. r22cr (a = L, R) describe spin de- 
pendent coupling of the dot QD2 to the ferromagnetic leads, 
whereas tj2cr is the hopping parameter between the j-th dot 
{j = 1,3) and the dot QD2. 



system of three coupled quantum dots, 



where, /^(e) 



^ \p{<^-t^c.)/kBT 



+ 1] ^ is the Fermi-Dirac 



distribution function in the lead a, G<(e) and G'"^''^(e) 
are the Fourier transforms of the lesser and retarded 
(advanced) Green functions of the dots, and r"(e) (for 
a = L,R) is a matrix which describes coupling of the 
dots to the leads. In the case under consideration the 
matrix r"(e) takes a simple form with only one nonzero 
matrix element. 








(4) 



Hqd eiadl^dia + 



J( = l,3) o- 



(1) 



where €ia is the energy level of the i-th dot {i = 1, 2, 3), 
and tj2(T is the hopping parameter between the j-th dot 



In the case under consideration one can write T22cr — 
r(l ± pl) for couphng of the dot to the left lead, and 
^22(7 — 7r(l ± Pr) (parallel configuration) or F^g^ = 
7F(1 ^ Pr) (antiparallel configuration) for coupling to 
the right lead. Here, tr =t (upper sign), a =| (lower 
sign), F is a constant, whereas 7 describes asymmetry in 
the coupling of the dot to the left and right leads. 

To calculate electric current one needs to know the 



{j = 1,3) and the dot QD2. Both e,^ and tj2a- are as- Green functions G 



r(a). 



sumed to be spin dependent in a general case. The sec- 
ond term of the Hamiltonian Hqd describes the intra-dot 
Coulomb interactions, with Ui (i = 1,2,3) denoting the 
corresponding Coulomb integrals. 

The last term of the system's Hamiltonian, Ht, de- 
scribes electron tunneling from the leads to the dot QD2 
(or vice versa), and takes the form 



2ka'-kQcr"2CT 



h.C. 



(2) 



(e). These can be found by the 
equation of motion method. First, we apply the equa- 
tion of motion to the causal Green function G'ycr(e). 
This equation of motion generates higher order Green 
functions, for which we write the corresponding equa- 
tions of motion. New higher order Green functions are 
then decoupled in terms of the Hartree-Fock decoupling 
scheme. This decoupling scheme, described in details in 
Refs^iii^, allows to express the higher-order Green func- 
tions through the lower-order ones. This, in turn, allows 
closing the relevant system of equations for the Green 
function Gija{e) and to write it in the matrix form as 



where V^ko- ^^'^ ^^'^ relevant matrix elements. Coupling of 
the dot to external leads can be parameterized in terms 
of r^(e) = 27: ^2l.'^2k;'5(e - Ckca). We assume that 
r"(e) is constant within the electron band, F"(e) = F" = 
const for e € {—D,D), and F"(e) = otherwise. Here, 
2D denotes the electron band width. Introducing the 
spin polarization pr of lead r {r — L,R) as Pr = {p^ — 
p~)/{p^ + p~), the coupling parameters can be expressed 
asr+^"^ = Tril±Pr), with F^ = {T+ + T-)/2. Here, p+ 
and p~ are the densities of states at the Fermi level for 
spin-majority and spin-minority electrons in the lead r, 
while F+ and F~ describe coupling of the dot QD2 to the 
lead r in the spin-majority and spin-minority channels, 
respectively. 
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where for simplicity the e-dependence of the Green func- 
tions and other parameters has not been indicated ex- 
plicitly (it will be restored where necessary). Apart from 
this, we have introduced the following notation: 



A,. 



n2is] 



Cjcr - Ui 



(6) 
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for 1 = 1,3, 



A2 = U2{e - e2a - U2) ^[tl2a{ni2S - n2is) 
+h2a{n^2S - "23?) + {Cls - C2s)n2S^a], (7) 



and 



B, = 1 



Utn2o 



U, 



(8) 



for i — 1,2,3. To simplify notation wc put here 

{nijcr) = {dfdj) = riija, {uic,) = (d^di) = n^cr, and 

C2a = J2ka^2ka{(iasd2s)- We have also defined the 
self-energy E^- as 



with 



(9) 



(10) 



for a = L,R. 

Having found the causal Green functions from equa- 
tion 5, one can calculate the retarded (advanced) Green 
functions as G[j^^(e) = Gijcr{e±iO~^)- The corresponding 
self-energy ^^'■''■'(e) = Eo-(e ± iO+) takes then the form, 



S;('^He)=A22.(e)T^r22., 



with 



and 



r22o- — ^22a + 



22CTJ 



A22.(e) = --r,2.1n 



(11) 



(12) 



(13) 



Now, we need to find the lesser Green function. This 
can be found from the corresponding equation of mo- 
tion, with the higher order Green functions calculated 
on taking into account the Langreth theorem^i and the 
Hartree-Fock decoupling scheme used when calculating 
the causal Green function. 



IV. NUMERICAL RESULTS 

Using the formulas derived above, we calculate numer- 
ically the basic transport characteristics, i.e. the conduc- 
tance and tunnel magnetoresistance. The latter quantity 
describes the resistance change when magnetic configura- 
tion of the system varies from antiparallel to parallel, and 
is described quantitatively by the ratio (i?AP — Rp)/ Rp, 
where i?AP and Rp denote the resistance in the antiparal- 
lel and parallel magnetic configurations, respectively. We 
assume that the dot levels are spin degenerate, = 
(for i = 1,2,3), and introduce two parameters, hi and 
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FIG. 2: Linear conductance in the parallel magnetic configu- 
ration, calculated as a function of the energy level 62 of QD2 
for indicated values of the intra-dot coupling parameters t2\ 
and t23 and for p = 0.4, = U^ = U^ = 0, hx = hj, = 0.2, 
7 = 1, and fcsT = 0.01. The energy is measured in the units 

of r. 



/i2, which describe separation of the levels t\ and 63 from 
the level £2, i.e. ci = £2 ^ ^1 and £3 = £2 + /i2- Moreover, 
we assume = = p and 7 = 1. 

Consider first the limit of vanishing Coulomb interac- 
tion, Ui — Q for i — 1,2, 3. In figure 2 we show the linear 
conductance in the parallel magnetic configuration, cal- 
culated as a function of the energy level £2 of QD2, mea- 
sured from the Fermi level of the leads. The total conduc- 
tance as well as that of individual spin channel are shown 
there. The energy level of the dot QDl is shifted down 
by 0.2 with respect to £2, while that of QD3 is shifted 
up by the same amount. Here and in the following the 
energy is measured in the units of F (F = F^ = T^). 
The hoping parameter t23 is constant while the param- 
eter t2i increases from zero in figure 2(a) to t2i = ^23 
in figure 2(f). Let us analyze now the main features of 
the conductance spectra. In figure 2(a) the dot QDl is 
entirely detached from the dot QD2 and therefore has no 
influence on transport. The other two dots are coupled 
and form bonding and anti-bonding like states. Due to a 
nonzero value of h2, these states are coupled differently 
to the leads. The linear conductance shows then two 
relatively broad peaks related to the bonding and anti- 
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bonding like states, with a deep minimum in between, 
where the conductance turns to zero due to destructive 
interference of electron waves transmitted via bonding 
and anti-bonding states. We note, that the broadening 
of peaks is determined mainly by the coupling to the 
leads, but temperature also plays a role. 

When the dot QDl becomes weakly coupled to the 
dot QD2, a narrow peak emerges in the conductance at 
the position corresponding to the energy level of the dot 
QDl (see figure 2(b-d)). The onset of this peak is a di- 
rect consequence of a nonzero value of ^21- The dots 
QDl and QD3 are then both coupled to the dot QD2 
and, disregarding coupling to the leads, one finds three 
new energy level of the whole three-dot system. The 
level corresponding to the new peak is only weakly cou- 
pled to the leads. When t2i increases further, this peak 
moves towards zero energy, and in the strictly symmet- 
rical case, t2i — ^23, the peak appears strictly at the 
position corresponding to the bare level of the dot QD2, 
£2 = 0. The spectrum becomes then symmetrical, which 
is a consequence of the symmetry in the system. The 
central peak, however, remains much narrower than the 
other two peaks. Similar scenario also holds for constant 
^21 and ^23 increasing from zero to ^23 = ^21- The new 
peak emerges now on the left side of the minimum and 
moves towards 62 = 0, as before. All the features dis- 
cussed above appear not only in the total conductance, 
but also in the conductance of each spin channel sep- 
arately, as clearly seen in figure 2. The heights of the 
side peaks in figure 2 are equal to 2 conductance quanta 
(e^//i), independently of the values of the inter-dot hop- 
ing parameters. The height of central narrow peak, how- 
ever, is smaller than that of the side broad peaks. 

The central narrow peak resembles the Dicke resonance 
in atomic physics, where a strong narrow emission line 
appears when the distance between atoms is smaller than 
the Fermi wavelength of the corresponding radiation. In 
our case the role of distance between atoms is plied by 
the distance in energy levels, hi and /12. When hi and 
/i2 increase, the width of the central line also increases. 
The two side peaks, in turn, become narrower and nar- 
rower. This is shown explicitly in figure 3, where the 
linear conductance for the symmetrical situation is pre- 
sented for increasing hi and /12 (the other parameters are 
the same as in figure 2). This behavior is quite reason- 
able, as the eigen-states of the dot system become more 
and more localized on the corresponding individual dots 
with increasing hi and /i2. Since the dots QDl and QD3 
are not coupled directly to the leads, the corresponding 
side peaks in the conductance become narrower and their 
width decreases with increasing hi and /i2- The central 
zero-energy level becomes then localized on the dot QD2, 
and the corresponding peak in the conductance becomes 
broader. Moreover, height of the central peaks increases 
and becomes equal to two conductance quanta, while that 
of the side peaks decreases with increasing hi and /i2- 

When hi = h2 ^ and ti2 = ^23 = ^ 7^ (as in 
figure 3), the two side peaks in the conductance appear 
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FIG. 3: Linear conductance in the parallel magnetic configu- 
ration, calculated as a function of the energy level 62 of QD2 
for indicated values of the level shifts hi and h2, and for 
t2i = t23 = 0.8, p = 0.4, = 1, Ui = U2 = U3 = 0, and 
ksT — 0.01. The energy is measured in the units of F. 



at the positions zt\/h^+21^. The third (central) peak 
occurs at e = 0, while the conductance reaches zero at 
±h. A particular situation occurs when hi — h2 — 0. 
The central peak, clearly seen in figure 3 for /ii, /12 7^ 0, 
does not occur then in the conductance. Since the system 
of three single-level QD's has three molecular-like states 
(denoted with the index 1, 2 and 3 for increasing energy), 
one could expect also three peaks in the conductance. 
However, the matrix elements of coupling between the 
molecular state |2) and the left and right leads vanish, 
the molecular state |2) becomes decoupled from the leads 
for hi = h2 = 0, which leads to the absence of the central 
peak in the conductance. As a consequence, one observes 
only two peaks in the linear conductance, situated at the 
positions e = ±V2t, and the transmission reaches zero 
at e = (eigenvalue of the state |2)). To show this more 
formally, let us perform the transformation of the dot's 
operators 



d2a 



V2 



{dia — d^a). 



(14) 



(15) 
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3a) 



(16) 



Assuming eia- — eo for i = 1, 2,3, the Hamiltonian of the 
isolated triple dots becomes then diagonal and takes the 
form 

Hqd = (eo + V2t)d\Ji„ + eodlj2a + (eo - V2t)dlJ^^. 

(17) 

In turn, the tunnelling Hamiltonian acquires then the 
form 



kacT i— 1,3 



where 



/2 



^3kcr — 2 2kcri 



and 



'^2kCT 



0. 



(18) 



(19) 



(20) 



(21) 



This clearly shows that the central molecular-like level 
becomes effectively decoupled from the leads. 

Let us consider now the influence of Coulomb cor- 
relation in the dots (nonzero parameters Ui for i — 
1,2,3) and nonzero temperature on the spectra described 
above. Figure 4(a) corresponds to figure 2(f), except the 
Coulomb parameters which in figure 4(a) are nonzero, 
Ui = U2 = = 1, while in figure 2(f) they vanish, 
Ui — U2 — U3 = 0. Owing to the nonzero Coulomb 
parameters, the number of peaks in the conductance be- 
comes doubled in comparison to that in the case with 
vanishing Coulomb interaction. The side peaks are broad 
as in figure 2(f), whereas all the other four peaks are 
very narrow. Figure 4(b) shows the conductance spec- 
trum for the same parameters as in figure 4(a), except 
temperature which now is higher. The peak heights are 
now smaller and the peaks become broader, as one might 
expect. Figures 4(c-e) show the conductance for a sys- 
tem with the Coulomb parameter vanishing for one or 
two dots. When the Coulomb parameter of a particu- 
lar dot vanishes, the number of peaks decreases by one. 
In all the cases, however, the side peaks are broad while 
the others are narrow. Figure 4(f) corresponds to fig- 
ure 4(a), but calculated for larger Coulomb correlation 
parameters, Ui — U2 ~ U3 — 4. The whole spectrum be- 
comes clearly split into two parts; one is similar to that 
for vanishing Coulomb parameters, and the second is its 
Coulomb counterpart, with a mirror plane somewhere in 
between the two parts. 

All the results presented above were calculated for the 
parallel configuration of the leads' magnetic moments. 
Qualitatively similar results have been obtained for the 
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FIG. 4: Linear conductance in the parallel magnetic configu- 
ration, calculated as a function of the energy level 62 of QD2 
for indicated values of temperature and intra-dot Coulomb 
repulsion parameters Ui,U2,Uz. The other parameters are: 
t2i = tn = 0.8, hi = h2 = 0.2, p = 0.4, and 7 = 1. The 
energy is measured in the units of F. 



antiparallel magnetic configuration. The quantitative 
difference leads to a nonzero TMR. In figure 5 we show 
TMR in the linear response regime as a function of the 
energy level 62 of QD2. When the dots' levels are well 
above or well below the Fermi level of the leads, the TMR 
tends to the JuUiere's value observed in planar magnetic 
tunnel junctions. When, in turn, the dot levels cross the 
Fermi energy, the situation becomes more complex and 
TMR significantly depends on the Coulomb parameters 
Ui. In the limit of C/j = for z = 1, 2, 3, see figure 5(a), 
TMR shows then a clear peak structure and fluctuates 
about the JuUiere's value. Qualitatively similar behavior 
of TMR was also found in transport through single quan- 
tum dot&'^^. Variation of TMR with the level position is 
more complex for nonzero values of the Coulomb corre- 
lation parameters, as shown in figure 5(b). Significant 
resonance-like peaks in the TMR occur now at certain 
values of the level positions, at which the conductance 
(current) is strongly suppressed by the interference ef- 
fects. Presence of such peaks in TMR is a consequence 
of the fact that the conditions for full destructive interfer- 
ence depend on magnetic configuration of the system for 
nonzero values of Ui. Outside the resonance-like peaks 
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FIG. 5: TMR in the linear response regime, calculated as a 
function of the energy level e2 of QD2 for indicated values of 
the intra-dot Coulomb repulsion parameters Ui,U2,U3, and 
for t2i = t23 = 0.8, hi = h2 = 0.2, p = 0.4, 7 = 1, and 
ksT — 0.01. The energy is measured in the units of F. 



the behavior of TMR is similar to that for vanishing Ui, 
except that the number of different peaks is larger due to 
the Coulomb counterparts in the conductance spectra. 



CONCLUSIONS 



In conclusion, we have calculated conductance of three 
coupled quantum dots. The conductance spectra, cal- 
culated in the Hartree-Fock approximation, clearly show 
presence of a central narrow peak, which resemblances 
the Dicke resonance in optical emission spectra of atoms. 
We have also shown, that the Coulomb correlations on 
the dots lead to doubling of the resonance peaks in the 
conductance. Additionally, we have calculated TMR 
effect associated with leads' magnetic moments rota- 
tion from antiparallel to parallel alignment, and shown 
that the intra-dot Coulomb correlation leads to some 
resonance-like enhancement of TMR at some values of 
the level positions. 
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